Time-dependent Ginzburg-Landau theory with floating nucleation kernel; 
FIR conductivity in the Abrikosov vortex lattice 
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We formulate the time-dependent Ginzburg-Landau theory, with the assumption of local equi- 
librium made in the reference frame floating with normal electrons. This theory with floating 
nucleation kernel is applied to the far infrared (FIR) conductivity in the Abrikosov vortex lattice. It 
yields better agreement with recent experimental data [PRB 79, 174525 (2009)] than the customary 
time-dependent Ginzburg-Landau theory. 
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The time-dependent Ginzburg-Landau (TDGL) equa- 
tion is a useful extension of the equilibrium Ginzburg- 
Landau theory. Unfortunately, microscopic derivations 
[J, [3, [3, 0, guarantee its validity under such restric- 
tive conditions that it seems more difficult to find justi- 
fied nontrivial applications than to solve it. The TDGL 
equation is thus most often applied beyond its nominal 
range of validity. 

As one leaves the familiar vicinity of the superconduct- 
ing phase transition and asymptotically slow processes, 
the intuitive foundation of the theory becomes shaky. 
The TDGL theory contains an assumption of local equi- 
librium, which is dependent on reference frame; when we 
adapt the equilibrium-based equation to non-equilibrium 
problems, we should at least work in the reference frame 
in which electrons are as close to local equilibrium as pos- 
sible. This is the frame floating with the normal current 
in the background of a superconducting condensate. To 
this end, in this paper we introduce what we refer to as 
a floating nucleation kernel. 

The standard TDGL theory is formulated using a ker- 
nel static in the laboratory system. We will show that 
compared to the TDGL theory in the floating system, the 
laboratory formulation lacks a term which is particularly 
important at high frequencies of the driving field. We will 
demonstrate the effects of this term on the conductivity 
in the sub-gap far-infrared (FIR) region. Comparing our 
results with recent FIR magneto-transmission measure- 
ments of Ikebe et al we will show that use of the 
floating nucleation kernel improves agreement between 
the theory and experimental data. 

Let us first describe the magneto-transmission mea- 
surement. It is performed on a thin layer perpendicularly 
penetrated by the magnetic field in the form of vortices. 
The incident FIR light is perpendicular to the surface 
and its electric field drives currents which determine the 
amplitude and phase of the transmitted light which is 
measured. 

Both the normal and the superconducting electrons are 
accelerated by the electric field and experience a fric- 



tion with the lattice. The friction of the condensate is 
much weaker since Joule heat develops only in vortex 
cores moving perpendicularly to the electric field. The 
relative contribution of these components to the current 
depends on the frequency of the driving field; the higher 
the frequency the higher will be the fraction of the nor- 
mal current. 

It is useful to inspect characteristic times for NbN, 
the material used by Ikebe et al The optical gap 
2 A = 5.3 meV implies the maximal sub-gap frequency 

< 10 THz. The mean time between two collisions of 
the normal electron is r„ ~ 5 fs, therefore during a single 
period of the sub-gap FIR field the electron loses momen- 
tum more than a hundred times. At zero magnetic field 
the condensate suffers no friction. The field of ampli- 
tude E accelerates the condensate to velocity e*E/ojm*, 
while a normal electron is accelerated to eETn/m. At the 
measurement temperature, T = 3 K and Tc — 15 K, the 
density of condensed electrons exceeds the normal den- 
sity, therefore the condensate clearly dominates the total 
current. A different situation obtains, however, for the 
Joule heat. The condensate current is out of phase with 
the driving electric field and generates no heat. The nor- 
mal current is in-phase, producing heat. If the magnetic 
field penetrates the sample, the condensate generates the 
Joule heat due to motion of vortices. We will see that 
for the sub-gap FIR frequencies the Joule heat value is 
much smaller than the amount of heat generated by nor- 
mal electrons. 

To identify the Joule heat, it is necessary to measure 
the transmission coefficient, including its phase. This al- 
lows one to determine the complex conductivity a with 
Im (7 giving the off-phase current and Re cr for the in- 
phase current. Ikebe et al 6] achieved this task by split- 
ting short pulses and mixing them again after one of 
branches passed through the sample. As mentioned, we 
will compare their experimentally established a with the- 
oretical predictions based on the TDGL theory in the 
laboratory and the floating coordinate system. 

We win use the electric field E(t) = Re [Ee"*'^^] and 
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current J(r) = Re [Je"*""^] . The complex conductivity 
is defined via J = crE. The current has a smaU Hall 
component which we neglect in our discussion for conve- 
nience. 

The TDGL equation derived using the static kernel ^ , 
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describes the evolution of the condensate including a re- 
laxation of the GL function towards its equilibrium 
value. The vector potential is that of the internal mag- 
netic field as well as the electric field of the FIR light ; 
B = V X A and E = — (1/c)9tA. The electric current 
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is composed of circulating diamagnetic currents and os- 
cillating response to the light. We solve Eq. ([1]) to linear 
order in E and eliminate the diamagnetic currents by av- 
eraging over the elementary cell of the Abrikosov vortex 
lattice; = (j^) = (B/^o) J^^^i^xdyis- The supercur- 
rent, = ct^E, gives the condensate conductivity 
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where t — T/Tc, b = B/Hc2 are the dimensionless tem- 
perature and magnetic field, cto is the normal state con- 
ductivity, /3a = 1.16 is the Abrikosov constant for hexag- 
onal vortex lattice, and Tg ~ r{l~t)/a. Deriving Eq. ([3]) 
we have used the GL parameter [Sl 
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The zero-temperature coherence length is determined 
by the upper critical field; = $o/(27ri7"2)- Here 
^c2 ~ 15 T is obtained via the linear extrapolation 
Hc2 = H^2i^ ~ from experimental data in Fig 3 of Q. 
The normal-state conductivity ctq = 2 • lO^/Jlcm, exper- 
imentally established at 20 K has weak temperature 
dependence and can be used at 3 K. 

In Fig. [T] one can see that the imaginary part of CTs from 
formula ([3]) reproduces recent experimental data of Ikebe 
et al [G;]. Here we use the GL parameter k — 40, the only 
fitting parameter in the present theory. It is adjusted to 
fit the imaginary part of the conductivity at 7 T. Our 
main interest is in the Joule heat given by the real part 
of the conductivity. 

Formula ([3]) was derived for the dense Abrikosov vor- 
tex lattice. Theoretically, the region of nominal validity 
is B > 4 T, at the temperature T = 3 K. It is there- 
fore somewhat surprising that theoretical curves of Im a 
slightly depart from the experimental data only at the 
lowest magnetic field i? = 1 T. 
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FIG. 1: Imaginary part of the conductivity giving non- 
dissipative currents: Thin lines are the superconducting con- 
densate conductivity Imas (dotted), the TDGL conductivity 
ImcTGL (full), and the two-fluid modification of the TDGL 
conductivity Imatf (dashed). The heavy line is the conduc- 
tivity Imo-fTj evaluated in the floating system. Experimental 
data of Ikebe et al [gj at 7 T (•) are in the nominal validity 
range of the TDGL theory, while the lower magnetic fields 
5 T, 3 T, and 1 T (o) are not. 



Due to the relaxation term Vdtip, the TDGL equa- 
tion ([1]) includes a damping and generates Joule heat 
y, Q = 4fcBTr(w/27r) (l^rV'P), where the brackets de- 
note the time average: (0) = (ci;/27r) Jg^'^^'^ dr^. The 
left-hand panel of Fig. [2] shows that the supercurrent 
produces Joule heat only at vortex cores. The right- 
hand panel of Fig. [5] presents the spatial distribution 
of the power absorbed by the condensate from the elec- 
tric field W — (js • E). The most intensive absorption is 
around vortices in regions elongated in the vertical direc- 
tion which is parallel to the electric field. Deep minima 
of the absorption are between vortices in horizontal rows. 
Comparing the two panels shows that the relation be- 
tween absorption and heat production is very non-local. 

The fraction of Joule heat due to the condensate is 
small. In Fig. [3] we compare the real part of the con- 
densate conductivity ([3]) with experiment. Indeed, the 
discrepancy between experimental data and RecTj indi- 
cates that the supercurrent produces only a minor part 
of the Joule heat; the normal current cannot be neglected 

From microscopic derivations 0, H S, of the GL 
theory it follows that the normal current and the super- 
current simply add. Adding the current J„ = a^iX + 
ir„u;)E which would appear in the normal state one ob- 
tains the TDGL conductivity 



ctgl = o-s 
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with the normal conductivity (T„ = cro(l -I- iTnUi). For 
experimentally established values ctq = 2 ■ lO^/Jlcm and 



3 




FIG. 2: Heat production (left) and the power absorption 
(right) in the hexagonal Abrikosov vortex lattice: Crosses 
denote centers of vortices. The electric field is polarised ver- 
tically so that vortices oscillate horizontally with amplitude 
shown by arrows. The Joule heat is produced at vortex cores, 
their horizontal motion is responsible for elongation of the 
heated region. Absorption of power is rather delocalised. Its 
maxima are also around vortex cores but elongated vertically. 
The rounded minima are between vortices. Difference of these 
two maps shows that the 'rigid' GL function transfers the 
power to be dissipated in cores. 
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FIG. 3: Real part of the conductivity giving Joule heat: 
Points are experimental data of Ikebe et al 6] for 7 T (•). The 
superconducting condensate contribution (dotted line) given 
by formula (|3]) is by an order of magnitude too small. The 
time-dependent Ginzburg-Landau theory (thin line) adds a 
contribution of normal electrons, see Eq. arriving at too 
high values. The two-fluid approach (dashed line) reduces 
the conductivity subtracting double-counted condensed elec- 
trons from the normal conductivity, see Eq. ((7)|. The float- 
ing kernel approach (heavy line) given by Eq. removes 
double-counting from the supercurrents and yields the closest 
agreement with experiment. 



T„ = 5 fs 6] , the normal conductivity yields a negligible 
contribution to ImcTGL, as seen in Fig.[l] but it provides 
the dominant contribution to RecrcL- One can see in 
Fig. [3] that Re ctgl is much closer to observed values than 
RecTs- It is higher than the observed values, however. 
This problem becomes more serious at lower magnetic 



fields, where the observed real part of total conductiv- 
ity is further reduced well below the level of the normal 
conductivity, see Fig. |31 while the TDGL conductivity is 
always larger, ImcrcL > Im(T„. 

The simple addition of normal current and supercur- 
rent works well close to the phase transition but it badly 
overestimates conductivity far from it. Apparently, it is 
insufficient simply to add the supercurrent and the nor- 
mal current; the electric field accelerates all electrons. 
Since electrons in the condensate escape frictional effects, 
this fraction of electrons must be removed in order to ob- 
tain the normal conductivity. An intuitive way to avoid 
double-counting of condensed electrons is to introduce 
a normal current reduced in the spirit of the two-fluid 
model, 
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The total current averaged over the elementary vortex 
lattice cell, J = + Jn, leads to a conductivity 



o-tf = as + {t + b) a„ 
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where we have evaluated the averaged normal fraction, 
1-2(|V'P) /n = t + b. One can see in Figs. [1] and [3] that 
the two-fluid conductivity yields the same non-dissipative 
currents described by Im atf as the TDGL theory, but 
that it allows for Retxtf smaller than the normal con- 
ductivity. In fact Re att is too small, when compared to 
experimental data. 

The reduced normal current ([6]) contradicts micro- 
scopic studies fl a a a i. 

Indeed, the total cur- 
rent is derived from the Nambu-Gor'kov Green function 
expanded in the gap, G ~ Gq + GqA*GqAGq, where 
Go gives j„ and the second term provides the supercur- 
rent. Apparently, the double-counting has to be remedied 
within the supercurrent itself. 

With this issue in mind we shift to our new formulation 
of the theory, expressing the nucleation of superconduc- 
tivity using the floating nucleation kernel. The Cooper 
pairs are created from electrons initially in the normal 
state, with mean velocity v = J„/(en). The free energy 
of condensation has to supply the kinetic energy which 
electrons gain going from the normal component into the 
condensate, therefore the stability condition reads 
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We note that quantum kinetic energy is in fact a non- 
local contribution of the nucleation kernel. For the float- 
ing kernel it depends exclusively on the velocity differ- 
ences of the normal and superconducting component 11 1 . 
The corresponding supercurrent 
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we can write as js 



e*v|(p|2 
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therefore this approach is free of double-counting. 

If an effect of velocity v on the GL function is neg- 
ligible, then if = ip and the total current jfk = js + Jn 
obtained with the floating kernel is not different from the 
current in the two- fluid approximation jtt = js + jn- In 
the presence of vortices, the kinetic energy is non-zero 
due to diamagnetic currents and the perturbation enters 
the TDGL equation in the linear order leading to changes 
of the GL function. The averaged total current Js + Jn 
then differs from J^-l- J„. The magneto-transmission thus 
allows us to test the TDGL theory formulated with the 
floating nucleation kernel. 

To obtain the conductivity we do not need to evaluate 
the modified GL function. The supercurrent modified 
by the inertial force m*dr'v is readily obtained from the 
condensate conductivity ([3]). The driving force in Eq. ([9]) 
is dr (— (e*/c)A — m*v) = e*E + z(tj/e*n)cr„E, therefore 
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The conductivity corresponding to the current 
given by 
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In Fig. [T] we compare Im (ifk with Im as ■ One can see 
that both values are very close except for at the small- 
est magnetic field where Im cTfk is closer to experimental 
data. 

In contrast, the Joule heat obtained within various 
approximations is rather different. In Fig. [5] we com- 
pare the standard TDGL theory with the floating kernel 
formulation. Although none of the approximations pro- 
vides satisfactory values, among the tested approaches 
our floating kernel prescription leads to values closest to 
experiment. 

In summary, we have formulated a version of TDGL 
theory using a floating nucleation kernel, meaning that 
the assumption of local equilibrium is applied to electrons 
in the moving reference frame of the normal current. 

When compared with standard TDGL theory in the 
context of far-infrared spectroscopy, we have found that 
the floating kernel formulation yields better agreement 
with experiment. In particular, recent published mea- 
surements of conductivity were considered; since we have 
established the GL parameter k from the non-dissipative 
response given by the imaginary part of the conductivity, 
our theory has no fitting parameters with respect to the 
Joule heat given by the real part of the conductivity. 

Finally, since use of this new approach does not gen- 
erally introduce significant additional complexity, it may 
be promising in the consideration of systems farther from 
equilibrium than is usually amenable to analysis via stan- 
dard TDGL theory. 
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FIG. 4: Real part of the conductivity giving the Joule heat: 
Points are experimental data of Ikebe et al |6|] for 7 T (•), 5 T 
(o), 3 T (□), and 1 T (A). The time-dependent Ginzburg- 
Landau theory (thin line) given by Eq. ([5} overestimates the 
dissipation. The two-fluid approach given by Eq. ([7} reduces 
the dissipation too much leading to the underestimate. The 
floating kernel approach ( heavy line) given by Eq. yields 
higher values although still smaller than experimental data. 
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